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SUMMARY 


OBJECTIVE 

To  illustrate  the  effect  of  update  randomization  on  the  envelope  statistics  of  K1  V- 
(if.N  reverberation  returns  and  to  examine  the  corresponding  probabilities  of  false  alarm  and 
detection  for  the  response  of  a broadband  energy  detector. 

RESULTS 

When  updates  are  limited  and  assigned  constant  magnitudes,  the  form  of  the  rever- 
beration envelope  is  altered  by  a reduction  in  the  probability  of  obtaining  high  intensity 
returns.  This  atfects  the  detection  statistics  for  a broadband  energy  detector  in  two  ways: 

1 . For  a given  signal-to-noise  ratio,  a signal  produced  with  updates  of  constant  mag- 
nitude will  result  in  erroneously  low  probability  of  false  alarm. 

2.  These  lowered  PFA  rates  cause  incorrect  placement  of  detection  thresholds 
which,  in  turn,  yield  erroneously  high  probability  of  detection. 


RECOMMENDATIONS 

Both  theoretical  calculations  and  computer  simulation  show  that  increasing  the  num- 
ber ot  updates  resolves  the  statistical  disparity.  However,  when  the  update  amplitudes  are 
randomized  according  to  a Rayleigh  distribution,  the  proper  statistics  are  obtained  independent 
of  the  number  of  updates  used  in  the  simulation. 
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INTRODUCTION 


The  rapid  development  of  microelectronic  technology  has  greatly  enhanced  the  signal 
processing  and  self-guidance  capabilities  of  lightweight  torpedoes.  To  accurately  evaluate  the 
performance  of  new  torpedo  designs,  including  such  features  as  matched  filter  processing  and 
coherent  signal  processing,  requires  a corresponding  increase  in  sophistication  ot  the  NOSC 
weapons  simulation  facility.  A major  element  of  that  improvement  is  the  realtime  implemen- 
tation of  REVGEN  (Reverberation  Generator),  a digital  simulation  of  ocean  reverberation. 


The  requirement  that  REVGEN  run  in  realtime  on  the  present  NOSC.  UNIVAC  1110 
based  hybrid  simulator  has  necessitated  several  simplifying  alterations  to  the  original  RI.V- 
GHN  concept.  This  note  examines  the  effects  of  these  simplifications  on  the  statistics  of  the 
REVGEN  output. 


The  physical  derivations  of  the  original  REVGEN  process  and  of  the  version  currently 

1 *>  } 

operative  on  the  NOSC  Hybrid  Simulator  have  already  received  ample  documentation.  — 
However,  for  the  purposes  of  a brief  review  and  to  establish  notation,  the  basics  of  this  process 
are  here  outlined. 


The  output  of  REVGEN  can  be  written  in  terms  of  the  following  matrix  equation. 

i.e.. 


Y = 

N-l 


Y 


■1 

L=0 


BJ.l.  XL 


(1  ) 


where  X represents  a sampled  vector  of  the  transmit  pulse  and  Y correspondingly  represents 
the  return  echo.  If  X and  Y sample  the  time  series  at  rates  equal  to  the  inverse  of  twice  the 
respective  band  widths  ( >v  y-,  \e  y)  of  the  transmit  pulse  and  return  echo,  it  can  be  shown  that 
the  B-matrix  contains  all  the  information  concerning  the  ocean  medium  that  is  available  to 
that  specific  ping-echo  process. ' 

The  /i-matrix  is  calculated  from  the  doppler-density  matrix  l)  in  the  following  wa 
N-l 

n - />  2njML/N 

H.!.L  ~ Z °J-L.Me 
\1=0 

where  V is  the  number  of  transmit  pulse  samples  in  X.  Lmax  = S - 1 . and  ./  indexes  the  time 
of  the  echo  return.  The  doppler-density  matrix  has  elements  that  correspond  to  different 
range  and  doppler  returns  in  the  echo.  The  rows  of  D j / represent  the  spectral  distribution 
of  the  scattercrs  at  a particular  range  slice  ./, 

To  conform  to  the  necessary  sample  rates  outlined  above,  the  adjacent  elements  of 
doppler-density  matrix  must  be  separated  by  1/2  w in  range  and  //_t/>  (t/>  = ping  length  tin 
frequency. 


C) 
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The  present  simulation  generates  reverberation  in  a three-step  process.  First  the 
doppler-density  matrix  is  calculated  from  data  defining  the  ocean  environment  and  the  tor- 
pedo state  (i.e.,  velocity,  beam  patterns,  etc.).  Second,  the  /(-matrix  is  formed  as  outlined  in 
Equation  (2).  Third,  the  echo  is  generated  from  the  /(-matrix  and  the  transmit  pulse  as  given 
in  Equation  ( 1 ). 

To  use  this  optimally  accurate  RFVGEN  calculation  in  the  weapon  simulator,  the 
rows  of  the  Dji  must  be  computed*  in  realtime.  For  typical  conditions  this  would  require 
256  rows  to  256  doppler-density  values  to  be  calculated  each  pulse  length  rp.  The  current 
capability  of  the  simulator  is  to  calculate  1 row  of  256  values  per  time  rp. 

The  disparity  between  the  optimum  number  range  slices  of  Dj^  per  pulse  length  and 
the  fact  that  only  one  can  be  calculated  has  led  to  certain  “update”  algorithms  which  serve 
to  fill  in  parts  of  the  doppler-density  matrix  without  resorting  to  the  lengthy  direct  calcula- 
tion of  the  ocean  algorithm. 

Figure  1 illustrates  how  this  update  procedure  is  accomplished.  Figure  1(a)  shows  a 
doppler-density  matrix  which  has  been  completely  filled  in  by  performing  the  ocean  algo- 
rithm calculation  for  each  range  slice  (separated  in  time  by  l/2w).  By  comparison.  Figure  1(b) 
shows  a similar  matrix  for  which  the  ocean  algorithm  has  been  called  only  once  per  pulse 
length. 

From  experience  it  is  known  that  the  average  magnitudes  of  each  element  in  the  rows 
of  Dji  do  not  change  much  during  one  pulse  length  (i.e.,  \Dj  ' \Dj+\  [ I where  .V  is  the 
number  of  range  slices  per  pulse  length). 

Since  the  average  magnitude  does  not  change  rapidly  with  range,  the  only  difference 
in  a particular  column  from  one  row  to  the  next  will  be  a randomization  in  phase  of  the  form 
ei<t>  where  <p  is  uniformly  distributed  between  zero  and  2ir.  Figure  1(c)  shows  a matrix  that 
has  been  filled  in  using  this  approximation.  The  first  row  in  each  pulse  length  has  been  calcu- 
lated using  the  ocean  algorithm.  The  rows  marked  with  a “u"  have  been  formed  using  the 
magnitudes  from  the  first  row  and  a random  phase  for  each  update  bin.  The  total  number  of 
rows  per  pulse  length  is  referred  to  as  the  number  of  updates  per  pulse  length.  The  current 
NOSC  realtime  simulator  is  capable  of  processing  8 updates  per  pulse  length.  Since  a ping 
may  typically  have  256  range  slices,  this  is  still  far  from  the  theoretical  optimum. 

*The  ocean  algorithm  which  is  used  to  construct  the  NOSC-REVGEN  doppler-density  matrix  partitions  up 
the  ocean  into  cells  corresponding  in  range  and  doppler  and  then  calculates  the  magnitude  ot  the  return 
from  each  cell  according  to  the  standard  sonar  equations. These  values  are  then  each  given  complex  ran- 
dom phase  and  loaded  into  the  appropriate  element  of  the  matrix. 

This  procedure  departs  from  the  true  form  of  the  doppler  density  in  that  it  assigns  a deterministic 
magnitude  (calculated  from  the  sonar  equations)  to  the  elements  of  Djp  while  in  reality  those  values  are 
stochastic  in  nature.  In  the  case  where  many  scatters  are  participating  in  an  instantaneous  echo  return  it 
has  been  shown  (in  Schwartz^  for  example)  that  the  magnitude  of  the  return  has  a Rayleigh  distribution. 

The  above  process  would  be  more  realistic  if,  after  the  mean  intensity  / //_  from  an  ocean  cell  was 
determined  from  the  sonar  equations,  the  magnitude  of  Djp  be  randomly  chosen  from  a population  of 
Rayleigh  random  numbers  having  the  mean  power  / y/. 
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B/N  < 1/2t 


(a)  THE  OPTIMUM  NUMBER 
OF  RANGE  SLICES 
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lb)  THE  NUMBER  OF  RANGE 
SLICES  RESTRICTED  BY 
DEMANDS  OF  REALTIME 
COMPUTATION 


c COMPUTED  VALUE 

u - UPDATE 

B - BANDWIDTH 

T * SAMPLE  TIME 
INTERVAL 


(c)  THE  INCREASED  NUMBER 
OF  RANGE  SLICES  DUE 
DUE  TO  UPDATE  APPROXI 
MATIONS 


figure  1.  Doppler-density  matrix  representations. 


At  this  point  several  important  questions  arise.  First,  how  will  this  under  sampling 
of  Dji  affect  the  results  of  a torpedo  simulation  run?  Second,  for  cases  where  the  change  in 
performance  is  significant,  is  there  a different  update  algorithm,  compatible  with  realtime 
operation  on  the  hybrid  simulator,  that  could  be  used  to  improve  performance? 

I • j 

The  answers  to  these  questions  are  quite  dependent  upon  the  type  of  signal  processing 

used  to  analyze  the  reverberation  return.  As  noted  by  Princehouse'  the  specifications  for  the 
bin  dimension  in  Dji  are  based  on  a matched  filter  processor  whose  transmit  pulse  simultane- 
ously contains  optimum  resolution  in  both  range  and  doppler.  In  situations  using  other 
processors,  a reduction  in  the  update  number  may  be  possible  without  affecting  the  output 
of  the  simulation  run. 

The  remainder  of  this  note  will  demonstrate  how  the  above  two  questions  can  be 
answered  by  studying  the  response  of  a simple  broadband  energy  detector  to  RHVGEN  out- 
put generated  using  several  different  update  algorithms.  It  will  first  be  shown  that  the  statis- 
tics of  the  REVGEN  signal-detector  process  can  be  represented  in  analytical  form.  The 
statistics  for  the  fully  updated,  properly  randomized  process  (i.e.,  real  world)  are  also  pre- 
sented in  closed  form.  From  these  expressions,  probability  of  false  alarm  statistics,  receiver 
operating  characteristic  curves,  and  power  curves  will  be  generated  to  demonstrate  the  effect 
of  the  update  algorithms  on  system  performance. 

These  results  will  also  be  demonstrated  using  actual  REVGEN  data. 


THEORY 


This  section  develops  techniques  to  statistically  characterize  1 ) the  reverberation 
return  from  point  scatterers  in  an  ideal  ocean;  2)  the  corresponding  REVGEN  output;  3)  the 
modified  REVGEN  output  aiu  its  application  to  the  theory  tor  two  dimensional  random 
walks:  and  4)  the  response  of  a broadband  energy  detector  to  real  and  simulated  reverbera- 
tion data. 

ENVELOPE  STATISTICS  FOR  A POINT  SCATTERING  MODEL 

The  physical  model  upon  which  REVGEN  is  based  assumes  that  the  reverberation  is 
the  sum  of  individual  echoes  from  point  scatterers  in  the  ocean  medium.  * The  expressions 
describing  the  reflected  return  from  such  a random  medium  of  point  scatterers  are  well 
known  and  adequately  reviewed  in  Ol'shevskiH  as  well  as  standard  texts  on  communications 
theory.  54' 

The  present  note  follows  chapter  5 of  Schwartz/’  to  introduce  the  statistics  of  signal 
propagation  through  so-called  random  or  fading  media. 


For  a medium  composed  of  randomly  placed,  stationary  point  scatterers  of  amplitude 
a i,  the  envelope  of  the  reverberation  signal  has  the  form 


- 1“L  <■* 


ml 


(3) 


where  U / is  determined  by  the  location  of  the  U1'  scatterer  in  the  medium.  Due  to  the  ran- 
dom placement  of  the  scatterers.  0 / is  a random  number  with  uniform  distribution.  When  no 
one  scatterer  dominates  the  return  ( i . e . . all  the  in  s are  of  comparable  strength ).  the  distribu- 
tion of  the  magnitude  of/-  can  be  shown  to  have  the  following  form. 


f/r) 


TL(,-/~/Jc 


(4) 


where 


8X0/  lias  the  following  properties 

1 . The  most  probable  value  of  r = a 

2.  /./ r)  = (ir/2)-  a 

3.  Median  value  of  /•=/./ $5  o 

•>  l 

4.  /./ r~ ) = power  = 2 a~ 


This  distribution,  named  after  Lord  Rayleigh  who  originally  solved  the  problem  in 


(5) 

<b) 
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5 . The  variance  E(  (E(r)  - rh ) = ( 2 - it  1 2 E 

6.  Probability  of  r exceeding  a threshold  rj  = e 


(9) 

(10) 


The  derivation  of  this  distribution  was  based  in  part  on  the  central  limit  theorem  and 
therefore  is  not  valid  for  cases  where  the  number  of  scatterers  is  small.  As  noted  by 
Princehouse,^  a typical*  sonar  might  ensonify  approximately  4 X 107  scatterers.  For  most 
pulse  lengths  this  would  insure  a large  number  of  scatterers  in  any  instantaneous  scattering 
volume. 

The  previous  paragraphs  describe  the  statistical  character  of  the  echo  return  from 
randomly  placed  scatterers  of  approximately  equal  amplitude.  Of  importance  to  detection 
theory  is  the  form  of  the  return  when  there  is  one  particular  dominant  component  in  the 
scattering  field  (e.g.,  a target). 

When  a component  of  amplitude  A is  added  to  the  previously  defined  group  of  scat- 
terers the  envelope  of  the  return  signal  has  the  following  distribution.9 


f(r) 


L_e-(r~  + A~  )/ 2 a~ 


!fj  < rA  /a~ ) 


(11 


where  Iq  is  the  modified  Bessel  function  of  the  first  kind  and  zero  order.  This  is  known  as  a 
Rician  distribution  and  can  be  seen  to  be  identical  to  the  Rayleigh  distribution  (Equation  (4)) 
when  A = 0. 


In  the  limit  that  rA  > o~  it  can  be  shown  that 
fir)  ^ e-<r-A)2/2o2 


(12) 


.’no* 


which  has  the  form  of  a Gaussian  distribution  centered  at  the  magnitude  of  the  specular 
component  A. 

Figure  2 shows  this  transition  from  Rayleigh  to  offset  Gaussian  distribution  by  plot- 
ting Rician  distributions  for  signal-to-noise  ratios  of  - 00  dB  (Rayleigh),  1 dB  (intermediate), 
and  10  dB  (off-centered  Gaussian).  The  signal-to-noise  ratio  (power)  has  the  form 


Power  target 
Power  Reverberation 


A2 


(13) 


The  foregoing  has  reviewed  the  statistical  nature  of  signals  that  would  be  received 
from  an  idealized  but  real  ocean  environment.  Next,  statistics  of  the  REVGEN  reverberation 
signal  will  be  determined. 


* 15-degree  beam  width  and  4000-metcr  range. 
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0.375 
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Figure  2.  Rayleigh  and  Rician  (off-set  Gaussian) distributions  for  SNR  of  1 and  10  dU. 

ENVELOPE  STATISTICS  FOR  KEVGEN 

Given  the  doppler-density  matrix.  REVGF.N  produces  a complex  time  series  (demod- 
ulated) corresponding  to  the  reverberation  echo.  From  Equations  ( 1 ) and  ( 2)  an  expression 
for  the  output  time  function  can  he  written  which  explicitly  shows  the  doppler-density 
matrix. 


G-I  2 '7.  v (14! 

1=0  M=0 

A few  observations  regarding  the  doppler-density  matrix  will  serve  to  simplify  the 
form  of  Equation  (12).  As  the  number  of  ocean  scatterers  falling  within  each  range  and  dop- 
pler  interval  becomes  large,  the  elements  /)/  / become  random  variables  of  the  form 

DjL‘djLJ*JL  ,15, 

where  <//  i is  distributed  according  to  a Rayleigh  distribution  and  0/  / is  a uniform  distribu- 
tion over  the  limits  10,  Jit).  This  follows  from  the  discussion  of  Rayleigh  statistics  presented 
previously.  For  the  case  of  a uniform  spatial  distribution  of  scatterers,  the  values  of/)/  / are 
statistically  independent. 


q 


L 


With  this  notation.  Equation  (14)  can  be  written 
X-l  N-l 


J=  X X xLdJ-LMe 


-J(  2wML/N  - 0 j_i  \f) 


( 1 6 ) 


L=0  M=0 


For  the  purpose  of  statistical  analysis  it  is  noted  that  the  exponential  factor  in  Equa- 
tion ( 16)  is  equivalent  to  e^L.M  where  U y ^ is  distributed  in  the  same  manner  as  <*>j  y. 


Making  the  definition 
aJ.L.M=xL  JJ-L.M' 

the  following  expression  may  be  written  for  Yj. 
y - - S n J°L-M 

L.M 

This  has  been  made  to  have  the  same  structure  as  Equation  (3). 


There  is  a fundamental  difference,  however,  between  the  values  aj  y j/  of  Equation 
( 18)  and  the  ay’  s of  Equation  (3).  The  subscript  L in  Equation  (3)  labels  each  scattering 
center  in  the  section  of  the  ocean  being  investigated  (as  noted  before,  this  is  a very  large  num- 
ber for  most  sonar  applications),  while  aj  y \j  labels  elements  in  a computer  simulation.  Due 
to  practical  considerations  such  as  the  time  and  cost  of  computation,  the  number  of  elements 
summed  in  Equation  ( 18)  must  be  kept  to  a minimum.  This  requires  that  one  ujl  \i  repre- 
sents the  scattering  from  many  of  the  ay's. 

If  when  computing  the  uj  y \f  s of  Equation  (17).  the  values  J j \j  are  chosen  accord- 
ing to  a Rayleigh  distribution,  then  |>yi  will  have  the  proper  statistical  distribution  indepen- 
dent of  the  number  of  elements  used  in  the  computation.  In  this  case.  Equation  ( 1 8)  is  the 
sum  of  complex  normal  random  variables  and  is  therefore  itself  a complex  random  variable 
having  magnitudes  distributed  according  to  a Rayleigh  distribution. 

However,  this  is  not  always  the  case.  When  the  elements  ay  y are  not  random  vari- 
ables but  are  fixed  according  to  a deterministic  calculation  (for  example  based  on  the  mean 
intensity  of  scatterers  from  within  a certain  volume),  then  the  only  randomness  in  >'./  occurs 
in  tfy  j\f.  This  in  turn  alters  the  statistical  distribution  of  \Yj\.  As  the  number  of  significant 
components  in  Equation  ( 18)  becomes  large,  \Yj\  will  approach  a Rayleigh  distribution. 
However,  when  the  number  of  components  is  small,  significant  differences  can  be  noted. 
From  Equation  ( 1 8 > it  can  be  seen  that  the  number  of  components  is  equal  to  the  product 
of  the  number  of  updates  per  pulse  length  and  the  number  of  non-zero  frequency  bins. 


Under  some  conditions  there  is  only  one  filled  frequency  bin  in  a reverberation  simu- 
lation. When  REVGEN  is  operating  in  its  update  mode  with  8 updates  per  pulse  length,  Yj 
would  be  formed  from  the  sum  of  8 complex  phasors  of  equal  magnitude;  i.e.. 
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(18a) 


o 

Yj=a  ^ e^‘ 

i=l 

where  0/  is  uniformly  distributed  on  |0,  lit).  In  this  situation  the  maximum  instantaneous 
power  ot  Yj  ( 6 4a- ) would  be  only  9 d B above  its  average  power  (Su-1. 

When  the  amplitudes c/y  / i/  are  fixed,  the  distribution  of  |)y|  is  equivalent  to  the 
distribution  ot  resultant  vector  length  in  a two-dimensional  random  walk  problem  where  a 
set  number  ot  stops  are  taken,  each  ot  a prescribed  length  a j [ \/  and  each  in  a uniformly 
random  direction  Oj  / . 

It  many  steps  ot  constant  length  are  taken  in  the  random  walk,  or  equivalently,  if 
many  updates  of  constant  amplitude  are  summed,  then  the  distribution  of  the  resultant  mag- 
nitudes approaches  a Rayleigh  distribution. 

Therefore  a Rayleigh  distribution  of  \ Yj\  can  be  achieved  in  either  of  two  ways: 

I ) by  using  many  updates  or  2)  by  randomizing  the  magnitudes  of  the  doppler-density 
matrix  elements  according  to  a Rayleigh  distribution.  (This  will  be  demonstrated  in  the 
results  section.) 

In  the  following  paragraphs  the  theory  of  the  two-dimensional  random  walk  will  he 
reviewed  to  provide  a basis  for  calculating  the  exact  statistical  nature  of  the  RKVGI  \ output 
when  the  amplitudes  ol  the  doppler-density  matrix  are  not  randomized. 


STATISTICS  OF  MODIFIED  REVGEN  (Two-Dimensional  Random  Walk) 

In  considering  the  REVGEN  statistics  ot  the  previous  section,  it  was  shown  that  mod- 
ification of  the  original  REVGEN  concept  also  alters  the  statistical  nature  of  the  resulting 
REVGEN  output.  In  effect,  the  distribution  of  the  magnitude  of  the  output  time  function 
I Y/\,  as  produced  by  the  modified  REVGEN  described  above,  is  equivalent  to  the  distribution 
of  the  resultant  vector  length  in  a two-dimensional  random  walk.  Thus  the  theory  to  be  pre- 
sented in  this  section  is  directly  applicable  to  REVGEN  output  when  the  amplitudes  of  the 
doppler-density  matrix  are  fixed  and  the  only  randomness  in  Yj  occurs  in  the  random  phase 
angles  0^  ^ 


Following  the  theory  presented  in  Rayleigh's  “Scientific  Papers."7  consider  the  proba- 
bility that  a traveler,  after  taking  //  steps  of  individual  lengths  u j.  a -> an.  in  directions  he 

chose  at  random,  will  consequently  arrive  within  a prescribed  distance  r from  his  starting 
point  (i.e..  Pn(  r.  <//,  a a,,)).  The  direction  of  his  first  st-p  is  plainly  a matter  of  indif- 

ference. On  (lie  other  hand,  at  his  next  n - 1 turns,  the  probability  that  the  angles  t)  lie 
within  the  limits  0 / and  0 / + JO /.  Oj  and  0 ->  + JO  ->....  d„_/  and  0n_j  + JOn_j  is 


— JO  i i 
Cv)"-1  ' 


10  -> 


JO 


II- 1 


(10) 


1 1 


The  probability  that  the  resultant  vector  is  less  than  r is  then  obtained  by  integrating 
the  above  probability  density  over  only  those  regions  ol  d space  (an  n-1  dimensional  niam- 
fold)  for  which 


1 > 

1 

;=  I 


i°i-i 


<r. 


The  evaluation  of  the  multi-dimensional  probability  integral 


Ue..Pn(r:a1.a2....an)-  J 

\ r — 71  ' ttuiu 


dO  j JO  i . . . dO 


niaa  <r 


n-i) 


(20) 


was  shown  by  Kluyver7  to  be  equivalent  to  the  following  definite  integral  in  one  variable. 
Kluyver’s  result  is 

OO 

Pn(r;  u j.  fli, . . . an)  = r J J jtrxIJ (fa jx)  J (fa2x*  ■ ■ ■ J(/anx * L,x 
0 


OO  )l 

= r f Jj(rx)  Jx  ]~[  (amx)  (21> 

0 m=l 

where  Jq  and  J / are  respectively  the  zero  and  first  order  Bessel  functions. 

Thus  the  problem  of  obtaining  numerical  values  tor  the  probability,  Pn,  lies  in  obtain- 
ing a method  for  accurate  numerical  evaluation  of  Kluyver’s  integral  (Equation  (21)).  The 
solution  of  this  problem  was  formulated  by  Bennett1 1 who  noted  that  for  the  special  case  ot 
n = i , the  problem  is  trivial;  for  the  case  of  n = 2 it  is  resolvable  in  terms  of  elementary  func- 
tions; and  for  the  case  of  n = 3 the  problem  can  be  evaluated  by  elliptical  integrals;  but  tor 
cases  of  n > 3,  an  alternate  method  of  evaluation  is  needed.  Bennett's  results,  summarized 
here,  conclude  that  the  probability  (Pn)  that  the  magnitude  of  the  resultant  vector  is  less  than 
r,  can  be  calculated  using  the  following  convergent  Fourier-Bessel  series  representation: 


Pn(r:al.a2,.  . . a„)  - 1 - j 


\1 


"m  c»i 


J(/rjm/A> 


m- 


I hn  ^ /h>^ 


(22) 


where  jm  is  the  mth  root,  in  ascending  order  of  magnitude,  of  the  equation  JQ(x)  - 0*  and 
where 


A = 


n 


1 

m=l 


m 


(23) 


See  Appendix  A for  table. 
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Or  conversely,  the  probability  that  the  magnitude  of  the  resultant  vector  will  "exceed’ 
r,  is  I minus  the  probability  given  in  Hquation  (22),  or  simply: 


CV<  i' s s 


■il 


I’m  1 ni  A t/'hn  * ^ 
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(2(i) 


I in  A / i 


Returning  to  the  case  of  the  traveling  pedestrian,  it  is  known  that  there  exists  an 
absolute  limit  on  the  distance  he  can  walk  in  n steps  since  each  has  a prescribed  length  um 
(the  limit  being  the  sum  of  the  lengths  of  all  n steps,  occurring  when  the  traveler  chooses  to 
walk  a straight  path).  Thus  it  is  impossible  tor  the  traveler  to  reach  a distance  greater  than 
the  sum 


I" 


in 


m=l 


Likewise,  in  hquation  ( 2b),  tor  the  case  where  r - . 1 . we  see  that  the  numerator  ot 
the  summing  term  includes  the  zero  order  Bessel  (unction  evaluated  at  jtn,  which  by  definition 
is  any  value  for  which  = <>.  Thus  when  r is  equal  to  the  limiting  value  A from  hqua- 

tion (23),  the  probability  that  the  resultant  vector  length  will  be  greater  than  r is  zero,  yield- 
ing precisely  the  result  expected  from  the  preceding  paragraph. 


From  these  results,  the  derived  probability  distribution  for  a finite  number  of  steps 
can  now  be  compared  to  the  Rayleigh  distribution  with  these  observations: 


(27) 


, y 

ver  = 2 a~  = 


2.  K(r~)  = power 

3.  With  target  component  of  amplitude  A , the  signal-to-noise  ratio  (power) 

*)  “> 

Power  target  A~  _ ,1~ 

' I ’ ' 


(28) 


Power  Reverberation 
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In  general,  this  derived  probability  distribution  lor  the  amplitudes  of  the  simulated 
reverberation  return,  generated  by  modified  REVGEN,  approximates  the  true  Rayleigh  dis- 
tribution witli  the  notable  exception  that  the  simulated  distribution  drops  to  zero  at  the 
point  which  represents  the  sum  of  the  component  magnitudes.  The  error  introduced  by  this 
"drop-off”  effect  is  essentially  negligible  for  large  numbers  of  components,  but  as  the  number 
of  significant  components  (i.e.,  updates)  is  reduced,  this  “drop-off”  error  becomes  increas- 
ingly significant.  It  is  this  discrepancy,  therefore,  and  its  resulting  effect  on  the  response  of 
a broadband  energy  detector,  that  becomes  the  subject  of  discussion  in  the  following  section. 


RESPONSE  OF  A BROADBAND  ENERGY  DETECTOR 

The  broadband  energy  detector  is  one  of  many  types  of  signal  processors  available  to 
current  sonar  systems  for  analyzing  reverberation  return.  This  simple  processor  examines  the 
envelope  of  the  reverberation  time  sequence,  testing  for  amplitudes  which  exceed  a predeter- 
mined "detection”  threshold.  Amplitudes  which  thus  exceed  the  set  threshold  are  sensed  and 
accordingly  recorded  as  target  detections.  The  following  section,  therefore,  examines  the 
statistics  of  this  "detection”  response  as  they  apply  to  the  probability  theory  already 
presented. 


From  the  first  part  of  the  theory  section  it  was  shown  that  for  an  ocean  medium 
where  many  seatterers  are  participating  in  an  instantaneous  echo  return,  the  magnitude  of  the 
return  follows  a Rayleigh  distribution^  and  when  an  additional  dominant  component  (e.g.. 
a target)  is  present  in  the  scattering  field,  the  envelope  of  the  return  then  follows  a Rician 
distribution.' ' Applying  these  results  to  the  broadband  energy  detector,  we  see  that  the 
probability  of  exceeding  a threshold,  T.  when  no  target  is  present  (thus  recording  a false 
alarm)  can  be  expressed  by  integrating  the  equation  derived  by  Rayleigh  from  T to  infinity. 
Thus  the  probability  of  false  alarm,  PFA.  for  a given  threshold  T is 


PFA  = 


(30) 


Likewise,  for  the  case  when  there  is  a target  present  in  the  scattering  field,  the  proba- 
bility that  the  magnitude  of  the  return  will  be  greater  than  the  threshold,  T,  (i.e.,  the  proba- 
bility of  detection,  PD)  is  expressed  by  integrating  the  Rician  function: 


OO 


In  general,  the  above  integral,  also  known  as  Marcum's  Q function,  has  presented 
some  difficulty  in  accurate  numerical  evaluation  for  non-zero  values  of  I . I bis  problem  and 
methods  for  its  solution  have  consequently  been  addressed  in  a variety  of  publications,  as 
summarized  in  Appendix  B. 
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At  tins  point,  Equations  ( 30)  and  (31)  provide  analy  tieal  expressions  lor  the  response 
statistics  of  the  broadband  energy  detector  when  it  is  operating  in  a fully  updated,  properly 
randonii/ed  system  (i.e.,  the  real  world).  Likewise,  the  third  theory  section  has  similarly 
presented  a well-defined  closed  form  expression  for  the  probability  of  exceeding  a threshold  . 
T.  that  applies  directly  to  the  simulated  REVGEN  return  when  it  is  operating  in  its  "update" 
mode.  Fq nation  (2b)  defines  the  probability  of  excess  applicable  to  both  the  condition  that 
no  target  is  present  (IM  A)  and  to  the  condition  that  there  is  a dominant  “target”  component 
present  in  the  scattering  field  (PD).  When  calculating  PFA,  the  sum  from  Kquation  (23) 
includes  only  the  appropriate  number  of  approximately  equal  scatterer  amplitudes,  but  when 
using  Kquation  ( 2b)  to  calculate  PI),  an  additional  dominant  amplitude  component  is  added 
to  the  Equation  (23)  sum. 

Although  Equation  ( 2b)  provides  a well-defined  representation  of  the  probability 
statistics  for  "update”  modified  REVGEN,  when  compared  to  the  probability  theory  of 
Equations  (30)  and  (31  ) for  a fully  updated,  properly  randomized  system,  apparent  differ- 
ences warrant  consideration  of  the  probability  “drop-off”  noted  earlier.  Some  questions 
remain,  such  as:  What  effect  do  the  limited  update  algorithms  used  to  generate  the  simulated 
REVGEN  output  have  on  the  corresponding  response  of  the  simple  broadband  energy  detec- 
tor? and  Mow  significant  is  this  effect  on  the  resulting  probabilities  of  false  alarm  and 
detection? 

Figures  3 and  4 represent  graphic  comparisons  of  the  “true"  Rayleigh  distribution 
and  the  simulated  distribution  of  REVGEN  output  generated  with  4 updates  per  pulse  length. 
Both  the  linear  ( Figure  3 1 and  the  logarithmic  ( Figure  4 ) graphs  show  evidence  of  the  proba- 
bility “drop-off"  previously  theorized.  When  concerned  with  obtaining  values  for  probabil- 
ities of  false  alarm,  we  are  specifically  concerned  with  the  area  in  the  tail  of  the  distribution 
beyond  the  value  of  the  detection  threshold  1 . I litis  it  is  clear  that  for  the  4-update  simu- 
lated REVGEN,  as  T approaches  the  drop-off  point,  the  values  for  PEA  drop  drastically  to 
zero. 

Similarly,  Figures  5 and  b illustrate  the  same  comparisons  made  in  Figures  3 and  4. 
except  that  the  simulated  RE  VGEN  return  was.  in  this  case,  generated  using  X updates  per 
pulse  length.  Linearly  (Figure  5),  this  simulated  distribution  appears  to  be  a much  more 
accurate  approximation  to  the  Rayleigh  distribution,  but  the  logarithmic  scale  of  Figure  b 
shows  that  although  the  drop  off  may  not  be  as  readily  apparent,  it  is  still  present  and  for 
higher  threshold  values,  it  will  still  affect  the  accuracy  of  probability  values  thus  obtained. 

In  general,  the  statistical  accuracy  of  simulated  probabilities  appears  to  be  greatly 
influenced  by  the  relative  proximity  of  the  detection  threshold  to  the  drop-off  point  for 
that  particular  simulation.  As  the  number  of  updates  used  to  generate  the  simulated  return 
increases,  the  value  of  the  drop-off  point  likewise  increases  and  the  error  (i.e.,  the  area  in  the 
tail  of  the  “true”  Rayleigh  distribution  beyond  this  point  I significantly  decreases. 

The  following  section,  therefore,  will  graphically  demonstrate  the  effect  of  various 
update  algorithms  on  probability  of  false  alarm  statistics,  signal-to-noise  ratio  power  curves 
and  receiver  operating  characteristic  curves. 
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RESULTS 


This  section  presents  graphic  illustrations  of  the  effect  of  various  update  algorithms 
on  the  probability  statistics  for  simulated  REVGFN  reverberation  return  and  the  correspond- 
ing detection  statistics  as  they  apply  to  the  various  outputs. 


PROBABILITY  OF  FALSE  ALARM 

The  probability  that  the  broadband  energy  detector  will  sense  a "target”  is  determined 
by  the  probability  that  the  amplitude  of  the  reverberation  envelope  exceeds  a set  detection 
threshold.  For  REVGEN  reverberation  returns,  the  update  algorithm  incorporated  in  the 
simulation  places  a maximum  on  the  possible  amplitude  of  the  return  envelope  and  because 
there  is  no  possibility  of  generating  amplitudes  greater  than  this  limit,  correspondingly  lowers 
the  probability  of  false  alarm  statistics.  Figure  7 shows  the  effect  various  update  algorithms 
have  on  their  respective  PFA  curves.  For  example,  at  the  PFA  level  of  10~3,  an  8-update 
REY'GEN  places  the  detection  threshold  at  approximately  5.8  signal-to-noise  power  (equiva- 
lent to  approximately  7.6  dB)  while  actual  Rayleigh  statistics  indicate  the  threshold  position 
of  7.0  signal-to-noise  power  (8.5  dB).  Considering  this  misplaced  threshold,  the  Rayleigh  sta- 
tistics show  that  this  7.6  dB  threshold  should  yield  an  actual  PF'A  closer  to  10“-  rather  than 
the  original,  erroneously  low  value  of  I0“3. 


(R- -2)  /2(SIGMA**2) 

I 1 1 T-1  ~ 1 1 I t 

-oo  6.02  9.03  10.8  12.0  13.0  13.8 

SNR  (dB) 

Figure  7.  PFA  versus  threshold  for  various  updates  per 
pulse  length;  also  PFA  for  Rayleigh  distribution. 
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In  addition  to  the  theoretical  curves  of  Figure  7,  corresponding  experimental  data 
(87,000  samples)  were  used  to  generate  Figures  8-15.  Applying  algorithms  for  2.  4,  and  8 
updates  per  pulse  length  with  fixed  amplitudes  and  random  phase,  actual  HFVCihN  output 
yielded  the  complex  time  series  of  the  reverberation  returns  shown  in  Figures  8-10.  For 
comparison.  Figure  I I presents  a reverberation  envelope  generated  with  both  random  phase 
and  Rayleigh-distributed  amplitudes.  Figures  I 2-15  then  use  this  generated  data  to  lend 
experimental  support  to  the  theoretic  probability  of  false  alarm  curves  from  Figure  7. 


With  a detection  threshold  set  at  3.5  dB,  Figures  8-10  illustrate  the  effect  various 
update  algorithms  (2,  4,  and  8 respectively)  have  on  corresponding  false  alarm  statistics.  I he 
apparent  amplitude  ceiling  shown  in  the  return  envelopes  of  Figures  8 and  0 results  in  unreal- 
istically low  PFA  rates  in  both  cases.  The  8-update  algorithm  of  Figure  1 0 shows  an  increase 
in  the  “detection”  peaks,  but  in  comparison  to  Figure  12  where  2 updates  per  pulse  length 
were  generated  with  amplitudes  randomized  according  to  a Rayleigh  distribution,  the  number 
of  false  alarms  for  8 updates  of  constant  amplitude  is  still  erroneously  low.  The  inaccuracies 
in  PFA  statistics  are  reiterated  in  Figures  1 2-14.  where  experimental  results  are  seen  to  agree 
with  the  previously  predicted  theoretical  PFA  values.  Of  particular  interest  are  the  experi- 
mental results  summarized  in  Figure  15.  A simulated  reverberation  return  with  update  ampli- 
tudes randomized  according  to  a Rayleigh  distribution  is  shown  to  produce  the  proper  PFA 
statistics  independent  of  update  number.  ( Figure  15  used  only  2 Rayleigh  randomized 
updates.) 


ENVELOPE,  dB  ENVELOPE,  dB 


Figure  8.  Experimental  representation  of  REVGEN  return  envelope  with  2 updates 
per  pulse  length,  constant  amplitudes,  and  uniform  phase.  Detection  threshold  set 
at  3.5  dB. 
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Figure  1 Experimental  PFA  results  (87,000  points  in  sample)  for  REVGEN  signal 
with  2 updates  per  pulse  length,  constant  amplitudes,  and  uniform  phase.  Solid  lines 
represen  theoretical  curves  for  2 updates  and  for  Rayleigh. 
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Figure  13.  Same  as  Figure  I 2 but  with  4 updates  per  pulse  length. 
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Figure  14.  Same  as  Figure  12  but  with  K updates  per  pulse  length. 
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Figure  I 5.  Same  as  Figure  I 2 except  the  update  algorithm  uses  2 updates 
per  pulse  length  having  Rayleigh  distributed  amplitudes. 
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RECEIVER  OPERATING  CHARACTERISTIC  CURVES 

Receiver  operating  characteristic  curves  graph  the  probability  of  false  alarm  versus 
the  probability  of  detection  for  various  detection  thresholds.  As  the  detection  threshold  is 
increased.  PE'A  decreases,  but  in  turn  PD  also  decreases.  ROC  curves,  therefore,  are  of  inter- 
est in  locating  that  detection  level  which  provides  a maximum  PD  for  a minimum  PEA.  Fig- 
ure 16  shows  the  Rayleigh  distribution  for  noise  alone,  the  Rician  distribution  with  a 10  d B 
target  signal  added,  and  the  areas  of  interest  in  determining  ROC  curves  < i.e. , shaded  PFA 
and  striped  PD). 

Given  the  altered  PFA  statistics  illustrated  in  the  previous  section,  it  is  clear  that 
RFVGEN  update  algorithms  will  likewise  alter  the  ROC  results  since  the  lowered  PFA  values 
are  used  in  the  ROC  comparison  of  PFA  and  PD.  Figure  1 7 presents  three  sets  of  ROC  curves 
(corresponding  to  signal-to-noise  ratios  of  8.  10,  and  12  d B » and  compares  those  values 
obtained  by  8 updates  per  pulse  length  simulated  REVC.EN  to  the  actual  ROC  values  from 
true  Rayleigh  and  Rician  distributions.  Because  the  REVGEN  PFA  drops  drastically.  Fig- 
ure 17  shows  that  for  small  increments  in  threshold  values  near  the  simulation’s  implied 
maximum,  PFA  decreases  rapidly  while  PD  changes  only  slightly. 

For  comparison.  Figure  18  presents  ROC  curves  similar  to  Figure  17  but  uses  a 
32-updates-per-pulse-length  algorithm.  In  all  three  cases,  the  increased  number  of  updates 
indicates  a closer  approximation  to  the  actual  Rayleigh/Rician  values. 

SIGNAL-TO-NOISE  RATIO 

This  third  type  of  graph  is  similar  to  the  ROC  curve  in  that  its  main  interest  lies  in  the 
relative  accuracy  of  values  obtained  for  the  probability  of  detection.  They  differ  in  that 
although  the  previous  ROC  curves  had  a fixed  signal-to-noise  ratio  (SNR)  and  a changing 
threshold,  the  following  plots  set  the  detection  threshold  at  a predetermined  level  and  vary 
the  signal-to-noise  ratio.  The  focus  of  Figures  19-21  is  to  indicate  that  at  various  SNR  levels, 
the  PD  obtained  with  various  versions  ot  REVGEN  simulations  are  erroneously  high  (the 
error  again  due  to  the  stochastic  nature  ot  the  update  algorithms  and  consequently,  to  the 
unreliable  values  for  PFA  which  cause  incorrect  placement  of  the  detection  threshold). 
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DISCUSSION 


From  the  preceding  section  it  is  seen  that  both  the  form  of  the  reverberation  signal 
and  the  related  false  alarm  and  detection  statistics  can  be  significantly  affected  by  changes 
in  the  update  algorithm. 

The  major  effect  is  the  reduction  in  the  probability  of  high  intensity  returns  in  the 
reverberation  signal  which  occurs  when  the  update  magnitudes  are  held  constant.  This  affects 
the  detection  statistics  in  two  ways.  First,  for  a given  SNR  threshold  the  probability  of  false 
alarm  will  be  less  for  a signal  produced  with  updates  of  constant  magnitude  than  for  the  cor- 
responding signal  with  updates  randomized  according  to  a Rayleigh  distribution. 

While  both  theoretical  and  experimental  evidence  indicate  that  increasing  the  number 
of  updates  resolves  this  disparity,  the  fact  remains  that  the  Rayleigh  randomized  updates 
produce  the  proper  statistics  (i.e.,  false  alarm  and  detection  statistics  for  a threshold  type 
energy  detector)  independent  of  the  update  number.  This  decoupling  of  the  detection  statis- 
tics from  the  update  algorithm  may,  in  some  instances,  permit  the  reduction  in  REVGHN 
computational  requirements  by  obviating  the  need  to  process  a large  number  of  updates  per 
pulse  length. 

While  it  has  been  shown  that  randomizing  the  update  magnitudes  decouples  the 
response  of  an  energy  detector  from  the  number  of  updates  used  to  form  the  REVGEN  signal. 

I it  should  be  concluded  that  in  tact  all  aspects  of  the  REVGEN  time-waveforms  become  inde- 

pendent of  update  number  when  the  proper  randomization  is  applied. 

Examples  of  this  occur  in  Figures  22-25  which  show  the  magnitude  of  the  base  band 
reverberation  (i.e.,  envelope  of  the  actual  signal)  for  two  pulse  shapes  and  for  two  different 
update  rates.  (In  each  figure,  5r  of  data  is  plotted.)  Figure  22  is  calculated  for  2 updates  per 
pulse  length  ot  Rayleigh  distribution.  Figure  23  represents  a signal  for  8 updates  per  pulse 
length  of  uniform  phase.  Both  figures  represent  the  return  from  a square  transmit  pulse.  For 
processors  which  classify  detections  based  on  a rapid  rise  in  signal  envelope  which  is  followed 
by  a sustained  high  value  for  some  period  of  time  (e.g..  t/2),  a decrease  in  the  number  of 
updates  might  trigger  more  false  alarms. 

Figures  24  and  25  (8  and  2 updates  respectively)  represent  the  reverberation  return 
when  a Chebyshev  window  has  been  applied  to  the  transmit  pulse.  Under  these  conditions 
the  envelopes  appear  more  realistic  than  for  the  square  pulse. 

The  conclusion  to  be  drawn  from  these  remarks  is  that  when  simplifications  are  made 
to  the  optimum  Rayleigh  simulation  (as  defined  by  the  pulse  length  and  band  width)  then  a 
validation  is  necessary  to  determine  the  response  of  the  specific  detection  logic  to  the  altered 
signal.  In  many  cases  the  validation  process  will  result  in  more  efficient  algorithms  for  gener- 
ating the  signal  while  at  the  same  time  insure  that  the  simplification  will  not  alter  system 
performance. 
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Figure  22.  Simulated  return  envelope  with  2 updates  per  pulse 
length,  uniform  phase.  Rayleigh  distributed  amplitudes,  and  a 
rectangular  transmit  pulse.  (Samples  taken  over  S pulse  lengths.) 
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Figure  23.  Same  as  Figure  22  with  8 updates  per  pulse  length  and 
constant  amplitudes. 
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Figure  24.  Same  as  Figure  22  (2  updates,  uniform  phase,  Rayleigh 
amplitudes)  except  with  a Chebyshev  window  applied  to  the 
transmit  pulse. 
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Figure  25.  Same  as  Figure  23  (8  updates,  uniform  phase,  constant 
amplitudes)  except  with  a Chebyshev  window  applied  to  the 
transmit  pulse. 
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Many  current  acoustic  detectors  are  based  on  a matched  filter  processor.  The  opera- 
tion ot  such  a tiller  is  very  ditlerent  Irom  that  of  a simple  energy  detector.  However,  it  can 
be  shown  that  the  statistical  problem  whic  h must  he  solved  in  order  to  evaluate  the  effect  of 
the  update  algorithm  is  ideiftieal  to  that  presented  in  the  theory  sec  tion  of  this  report. 

I he  details  ol  the  application  ol  this  theory  to  the  problem  of  a matched  filter  proc- 
essor and  the  ettect  ol  various  update  algorithms  is  the  topic  of  a future  report. 
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APPENDIX  A.  ZEROS  OF  THE  BESSEL  FUNCTION  J 


s J()'s 

I 2.404825557605772768621 63 1 879 
3 8.65372791 201 101221695419871  I 
5 1 4.9309 1770848778594776259399 
7 21.21 163662987925895907839339 
9 27.49347913204025479587728828 

1 1 33.77582021357356868423854639 
13  40.05842576462823929479930742 
15  46.341 18837166181401868578892 
17  52.6240518411 1499602925128542 
1 9 58.906983926080942 1 3283440668 

2 I 65.1 8996 480020686044063603378 
23  7 1 .47298 1 60359373282506307389 
25  77.75602563038805503773937193 
27  84.039090776938 1 90 1 5787963838 
29  90.322 1 726372 104800557 1 776680 
3 1 96.605267950996268778 12161 736 
33  102.8883742541947945964200347 
35  1 09. 1 7 1 4896498053835520659775 
37  115.4546126536669396281 177571 
39  1 2 1 .7377420879509629652343629 
41  128.0208770060083240797635598 
43  134.3040166383054660993528973 
45  140.5871603528542965484889028 
47  146.8703076257966495941327053 
49  153.1534580192278924875916260 
5 1 159.43661 1 1 64263 1 4632349 1 0383 
53  165.7197667479550208666904272 
55  1 72.0029245030782002 1 540043 1 9 
57  178.2860842000737706814839741 
59  184.569245640638718141 1 154117 
61  190.8524086525815223217807994 
63  1 97. 1 3557308566 1 4 1 473621 2003 1 
65  203.4 1 8738808 1 98646 1 7 1 2488207 
67  209.7019057042940751974824471 
69  2 1 5.98507367 1 5340 1 3 1 569956383 
71  222.2682426190843143412845248 
73  228.551412466098813301 197841 1 
75  234.8345831403832410198080237 
77  241.1 177545772680225149969715 
79  247.400926718652824848087941  1 
8 1 253.6840995 1 2 1 9308 1 0046928089 
83  259.9672729106044715730647780 
85  266.2504468710658801  188409547 
87  272.53362 135470493 1 4535  1 43489 
89  278.8 1 6796326 1 5308657846 1 7743 
9 1 285.09997 1 753 1 59564539 1 305882 
93  291 .383 14760625521 22416890505 
95  297.6663238584589425240491906 
97  303.949500485020581  1060577802 
99  310.2326774631949609525960935 


s Jo's 

2 5.5200781 102863106495966041 12 

4 1 1.79153443901428161374304491 

6 1 8.07 1 0639679 1 092254314788297 

8 24.35247 1 5307493027370579448 1 

10  30.6346064684319751 1754957897 

1 2 36.9 1 709835366404397976949305 

1 4 43. 19979 17131 7673035752407277 

1 6 49.4826098973978 1 7 1 7360276 1 58 

IX  55.76551075501997931 168349281 
20  62.04846919022716988285250031 
22  68.33 146932985679827099230388 

24  74.61450064370183788382054050 

26  80.897555871 13762786377214353 

28  87.180629843641 15365126180505 

30  93.4637 1 878 1 94477417 1 1 9059 1 58 

32  99.7468 1 985868059647027997904 

34  1 06.0299309 164516155101 769 1 76 

36  1 1 2.3 1 30502804949096274945066 

38  1 18.5961766308725317156293849 

40  1 24.8793089132329460452591288 

42  131.1624462752139146078961 154 

44  137.4455880202842777X77826619 

46  143.7287335736897325339507184 

48  150.01 18824569547574908804734 

50  156.2950342685335238195494957 

52  1 62.578 1 886689466775 1 90598094 

54  168.8613453692358256874563379 

56  1 75.1445041 2 1 9027430653727983 

58  181 .4276647 1 373 1 0507942 1 22000 

60  187.7108269600493597800752662 

62  193.9939907001091 197899403972 

64  200.27715579333241 17833620879 

66  206.5603221 162444736554572454 

68  2 1 2. 84348955  994948  2 7 5074  5 5 66 1 

70  219.1 26658028040567465 1891314 

72  225.4098274348593298985131824 

74  23 1 .6929977040385387809737966 

76  237.9761687672756628555575638 
78  244.2593405632956825588065766 

80  250.5425 1 30369699554704463657 

82  256.8256861385644130243177332 

84  263.1088598230954706927040030 

86  269.392034049776067 1 38 1 733966 

88  275 .67  5 208  78 1 537453847897  7954 

90  28 1 .9583839846 149 1 98543 167101 

92  288.2415596281876964381 173954 

94  294.5247356840649514582331775 

96  300.8079 1 2 1 264 1 1 1 3477 1 6689693 

98  307.0910889315050391 147737974 

100  313.3742660775278447196902455 
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s Jo’s 


s Jo’s 


01  316.5158547720429222184775231 
03  3 2 2 .79903  2 39  2 25 5 5 5 2003  247  545  3 
05  329.08221 0305998 5 5 80383430203 
07  335.3653884967741359209434008 
09  341.6485669492980779584522592 
1 1 347.9317456493901970068081779 
13  354.2149245838764120723287290 
1 5 360.498 1 03740501 07  295 1 4569438 
17  366.781283107848296151 1060446 
1 9 373.06446267527 1 2497490 1 83464 
2 1 379.3476424328284656 1 30799003 
23  385.6308223712263775181 195946 
25  391.91400248176738648721 14046 
27  398.1971827563028428764478006 
29  404.480363 1 87 1 904105848777601 
31  41 0.7635437672553441 9064 1 5374 
33  417.0467244897552663636082517 
35  423.32990534834808 1 9004966430 
37  429.6 1 30863370627072752 1 92287 
39  435.8962674502723316245713017 
41  442. 1 7944868266995738 16391445 
43  448.4626300292459969898 1 67  807 
45  454.74581 14852677268416657443 
47  461.0289930462604212672619595 
49  467.3 1 2 1 747079900084536543870 
5 1 473.595356466447 1 069586500785 
53  479.8785383178323162862533058 
55  486.1617202585426480732848010 
57  492.4449022851589960166860096 
59  498.7280843944345529380228240 
61  505.01 12665832840924990941259 
63  51 1.294448848774041 1915926959 
65  517.577631 1881 132734468657438 
67  523.8608135986445691552365105 
69  530. 14399607783667864048 1 7757 
71  536.42717862327694528441 19830 
73  542.7 1 036 1 2326644406093 1 89 1 59 
75  548.9935439038035707637447593 
77  555.2767266345981 170733475001 
79  561.5599094230456766607358788 
81  567.8430922672324721474293356 
83  574. 1 26275 1 65328502 1 639868063 
85  580.4094581 155830068429281626 
87  586.692641 1163202246816488596 
89  592.9758241659354191641086158 
91  599.259007262891 1553427808792 
93  605.5421904057138082257789970 
95  61  1.82537359299028630561 10142 
97  618.1085568233649549210978179 
99  624.39 1 7400955367453763497779 


102  31 9.657443544376 1 599488873476 
104  325.9406213134966851673704894 
1 06  332.2337993677396373454899 187 
108  338.5069776912284979209794860 
110  344.7901562692440017789753615 
112  351.0733350881205833077719639 
1 14  357.3565141351537464179662541 
1 1 6 363.6396933985 1 7037323 1 1 1 9046 
118  369.9228728671874783981774360 
1 20  376.206052530878473967667 1 365 
122  382.4892323799793288604465262 
124  388.7724124055006316474345647 
126  395.0555925990248496581620643 
128  401.33877295266156464175501 10 
1 30  407.62 1 953459006848364 1 8 1 593 1 
132  413.90513411 1 1063382384651423 
1 34  420. 1 8831490242 1 62570666570 1 8 
136  426.4714958267996157372832039 
138  432.7546768784445554823246178 
140  439.0378580518924647081334559 
142  445.3210393419877307942659586 
144  451 .60422074386165749592797 1 7 
146  457.8874022529127798062941462 
148  464.1705838647887775767442078 
1 50  470.453765575369838432 1 633945 
1 52  476.7369473807533362768792329 
154  483.0201292772397056003975080 
156  489.30331 12613194040977457813 
158  495.5864933296608670207868442 
160  501 .869675479099366350 1078187 
1 62  508. 1528577066266964738964728 
164  51 4.43604000938 161571 20541 695 
1 66  520.7 1 92223846409798446849400 
168  527.00240482981 15098931093407 
170  533.2855873424221418449084152 
172  539.5687699201 169108873902364 
174  5 45. 8519  5 25  6064832708195  3 1724 
176  5 52. 1 35 1 3526 1 87 1 2033325405405 
178  558.4 1 83 1 802 1 7369000250756037 
180  564.7015008382879538856080555 
182  570.9846837096530614617208419 
184  577.2678666340423902087180161 
186  583.551049609743192489371 1029 
188  589.8342326351 156999001434814 
190  596.1 174157085892772413556898 
1 92  602.4005988286588 1 7 1 745 1 75769 
1 94  608.68378 1 99388 1 358 1 754902234 
1 96  6 1 4.9669652028729098 1 363 1 9040 
198  621.2501484543054706791861 193 
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APPENDIX  B 


In  radar-detection  theory  and  in  other  problems  involving  carrier  signals,  it  is  often 
necessary  to  compute  Marcum's  Q Function.  The  Q function  is  defined  as: 


(B'l ) 


or  more  simply 


<J(a,b)=  J v exp  j /(/ax)d.v  ( B-2 ) 

/> 

where  !q  = zero  order  Bessel  function. 

While  there  do  exist  other  asymptotic*  and  recursive-  formulas  for  solving  the  0 
function,  the  recursive  method  derived  by  McGee appears  to  be  more  accurate  and  generally 
the  most  stable. 


This  recursive  method  can  be  derived  by  replacing  the  zero-order  Bessel  function  by 
its  series  expansion. 


_/ 

(>i!)2 


(B-3) 


inverting  the  order  of  summation  and  integration,  and  then  rearranging  the  Brennen  and 
Reed  summation-  to  be  a power  series  in  b instead  of  in  a.  Thus  we  obtain: 


(B-41 


(tK-J 

y* = ~rrr~ 


hk=,ik~i  +Jk 


d (f)^ 

d*=— ~k — 


After  N iterations,  the  error  in  Q is  given  by  the  remainder  term: 


E=  J fk,lk 


To  determine  a bound  on  E,  we  consider  the  iterative  equations  above,  which  when 
algebraically  substituted  within  one  another,  result  in: 


fk+1  ,lk+l 


\b2(2k  + a2i\ 
L 4k(k+2)  J J 


thus  giving  the  ratio  of  succeeding  terms  as 


fk+1  ,lk+l 


[b2<  2k  + a2) 
L 4k(k  + 2)  . 


(k  > l ) 


Note  that  as  k increases,  this  ratio  decreases.  Therefore  recursive  testing  for  the  con- 
dition that  this  ratio  of  successive  terms 

(fk+1  hk+l\  _ 1 

v-&rn 

will  then  imply  that  E < fxbx 

Thus  the  summation  can  be  terminated  when  two  conditions  are  met: 

1 . .he  ratio  of  succeeding  terms  < 1/2  (implying  E < f h ) 

2.  fx^x  *s  determined  to  be  an  insignificant  subtraction  from  the  existing  value  for  Q. 


36 


REFERENCES 


B-l.  ('.  W.  Helstrom.  “Statistical  Theory  of  Signal  Detection.”  2nd  eel . New  York  Pergamon. 
1 l)(i8,  p.  45  I . eq.  F-l  6. 

B-2.  L.  T.  Brennan  and  1.  S.  Reed,  "A  Recursive  Method  ot  Computing  the  (,)  Function,” 

IEEE  Trans.  Information  Theory  (Correspondence),  vol  I I -I  I . April  I ‘ > « > 5 
pg.  312-313. 

B-3.  W.  F.  McGee,  "Another  Recursive  Method  of  Computing  the  (,)  I unction  II  I I I rails 


